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^ ■ Abstract 

We solve the quantum Vlasov equation for fermions and bosons, incorporating 
^ ■ spontaneous pair creation in the presence of back-reactions and colhsions. Pair 

, creation is initiated by an external impulse field and the source term is non- 

i Markovian. A simultaneous solution of Maxwell's equation in the presence 

I of feedback yields an internal current and electric field that exhibit plasma 

oscillations with a period Tpi. Allowing for collisions, these oscillations are 
■ damped on a time-scale, Tr, determined by the collision frequency. Plasma 

oscillations cannot affect the early stages of the formation of a quark-gluon 
plasma unless Tr ^ Tpi and Tpi ~ 1/Aqcd ~ Ifm/c. 

O 

Ch ■ I. INTRODUCTION 



^ . The nonperturbative Schwinger mechanism which describes the spontaneous forma- 

tion of fermion-antifermion pairs, has been used 0] to model the formation of a quark-gluon 
plasma in heavy ion collisions. In this approach nucleon-nucleon collisions lead to the cre- 
ation of fiux-tubes, in which quark-antiquark pairs are connected by a strong colour-electric 
field. The energy density (string tension) acts like a strong background field and particle- 
antiparticle pairs are created via the Schwinger mechanism These charged particles 

polarise the vacuum and are accelerated in the external field. Their motion generates a field 
that in turn modifies the initial background field and, in the absence of further interactions, 
that back-reaction induces plasma oscillations. 

The back-reaction phenomenon has become a focus of attention in recent years, both in 
general and as it can arise in the pre-equilibrium stage of an heavy ion collision. Theoretical 
approaches as diverse as field theory P-[TO[| and transport equations |1TT|JT^ have been applied. 



The link between treatments based on the field equations and the formulation of a Boltzmann 



equation was recently investigated |T^,|TJ]. These studies show that the resulting kinetic 
equation has a non-Markovian source term. For weak fields there is no overlap between the 
time-scales characterising vacuum tunnelling and the period between pair production events: 



Tqu, Tprod, ciud the Markoviau approximation to the quantum Vlasov equation is valid []14 
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However, for strong background fields there is an overlap between these time-scales and this 



makes the non-Markovian nature of the source term very important |115HT7[| . Back-reactions 
and collisions introduce at least two more time-scales: the plasma oscillation period, Tpi, and 
the collision period, r^, and their impact is an integral focus of this article. Furthermore, in 
contrast to other recent studies we induce particle production by a time- dependent 



external field. 

In Sect. II we review the main equations and results for particle creation using a non- 
Markovian source term. In Sect. Ill we derive the renormalised Maxwell equation determined 
by the external and internal fields and, for special choices of the external field, present 
numerical results obtained by solving the coupled system of kinetic and Maxwell equations 
for bosons and fermions, with and without a simple collision term. We summarise our results 
in Sect. IV. 



II. PAIR CREATION WITH A NON-MARKOVIAN SOURCE TERM 

We consider an external, spatially-homogeneous, time-dependent vector potential A^, in 
Coulomb gauge: Aq = 0, and write A = (0, 0, A(t)). The corresponding electric field is 

E{t) = -A{t) := -dA{t) /dt . (1) 

The kinetic equation satisfied by the single-particle distribution function: f± ("+" for 
bosons, "— " for fermions) is 

^^ = S^{p,t), (2) 
where the source term is momentum- and time- dependent: 

S±ip, t) = lw±{t) f dt'Wi{t')F^{p, t) cos[x(t', t)] , (3) 

2, J —oo 

with x(t',t) := 2[9(t) — 6(t')] describing the difference between the dynamical phases 

Q{t) = f dt'uj{t') . (4) 

J —oo 

Here the total energy is 



u;{t) = ^el + Plit) (5) 

where e±_ = \J m? + is the transverse energy and we have introduced the kinetic momen- 
tum: P = {p±,Pii{t)), with p± = (pi,p2), P\\{t) = p\\ - eA{t). 

Equation (0) was recently derived from the underlying quantum field theory |T3Hl5| and 



exhibits a number of interesting new features. For example, in a strong background field 
its solutions describe an enhancement in the boson production rate and a suppression of 
fermion production. There are two aspects of Eq. (H) that generate such differences 



between the solutions for fermions and bosons: the different transition coefficients 
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where the degeneracy factor is (?+ = 1 for bosons and = 2 for fermions; and the statistical 
factor: F±{p,t) = [1 ± 2/±(p, t)]. 

The kinetic equation, Eq. (^, is non-Markovian for two reasons: (i) the source term on 
the right-hand-side (r.h.s) requires knowledge of the entire history of the evolution of the 
distribution function, from t^oo ~^ t; and (ii), even in the low density limit {F{t) = 1), 
the integrand is a nonlocal function of time as is apparent in the coherent phase oscillation 
term: cos[x(t', t)]. The mean field approaches of Refs. [p|-|TO[] also incorporate non-Markovian 
effects in particle production. However, the merit of a kinetic formulation lies in the ability 
to make a simple and direct connection with widely used approximations. 

In the low density limit the source term is independent of the distribution function 

Slip, t) = ^W±(t) r dt'W±{t') cos[x{t', t)] (7) 

Z J —oo 

and Eq. (g) becomes 

mn^suM. (8) 

(The low-density limit can only be self-consistent for weak fields.) Even in this case there are 
differences between the solutions for fermions and bosons because of the different coefficients 
W±(t), and the equation remains nonlocal in time. Equation (|^) has the general solution 

fl{p,t)= f dt'Sl{p,t'), (9) 



which provides an excellent approximation to the solution of the complete equation when 
the background field strength is small compared to the transverse energy. 

The ideal Markov limit was found in Ref. fl^, where a further asymptotic expansion 
was employed and a local source term for weak electric fields was derived. In this case 
Tqu < Tprod- Howcvcr, for Very strong fields a clear separation of these time-scales is not 
possible and the kinetic equation must be solved in its non-Markovian form where memory 
effects are important [|1^ . 



Equation (^ is an integro-differential equation. It can be re-expressed by introducing 
v±{p, t) = f dt' W±(p, t')F^{p, t') cos[x(p, t, t')] , (10) 
z±ip, t) = f dt' >V±(p, t')F±{p, t') sin[x(p, t, t')] , (11) 



in which case we have 



^^^ + ei^(t)^^^ = iW,(P,t').,(P,t), (12) 
+ eE{t) ^^^^^ = W±(P,t')i^±(P,t) - 2u;iP)z^iP,t) , (13) 
^^ + eE(t)^^^2.(PK(P,t), (14) 
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with the initial conditions /±(to) = v±(to) = -2±(^o) = 0, where to — oo. This coupled 
system of linear differential equations is much simpler to solve numerically. 



III. BACK-REACTIONS 

A. The MEixwell equation 

In recent years the effect of back-reactions in inflationary cosmology and also in the 
evolution of a quark-gluon plasma has been studied extensively. In both cases the particles 
produced by the strong background field modify that field: in cosmology it is the time- 
dependent gravitational field, which couples via the masses, and in a quark-gluon plasma, 
it is the chromoelectric field affected by the partons' colour charge. 



In our previous studies we have considered constant [|T^ and simply-constructed time- 
dependent [0 Abelian electric fields but ignored the effect of back- reactions; i.e., that 
the particles produced by the background field are accelerated by that field, generating a 
current that opposes and weakens it, and can also lead to plasma oscillations. The effect 
of this back-reaction on the induced field is accounted for by solving Maxwell's equation: 
E(t) = —j{t). Herein we assume that the plasma is initially produced by an external field, 
Eexit), excited by an external current, jexit), such as might represent a heavy ion collision 
and this is our model-dependent input. The total field is the sum of that external field 
and an internal field, Ein(t), generated by the internal current, jmit), that characterises the 
behaviour of the particles produced. Hence the total field and the total current are given by 

E{t) = EUt) + Eex{t) , (15) 
Jit) =jUt)+ Jexit) . (16) 

Continued spontaneous production of charged particle pairs creates a polarisation cur- 
rent, jpoiit), that depends on the particle production rate, Sip,t). Meanwhile the motion 
of the existing particles in the plasma generates a conduction current, jcondit), that de- 
pends on their momentum distribution, /(p, t). The internal current is the sum of these two 
contributions 

Einit) = -jin = -jcondit) - jpoiit) ■ (17) 

At mean field level the currents can be obtained directly from the constraint of local energy 
density conservation: e = 0, where 

eit) = lE\t) + 2 J ^,^iP, t)fip, t) . (18) 



For bosons this constraint yields 



(27r)3 ujip,t) 
and we can identify the conduction current 



^^iP,t) dfip,t) 



ujiPjt) dt 



4 



(20) 



and, using Eq. (|), the polarisation current 



jpol{t) 



E{t) J (27r)3 

Thus, using Eqs. (|TUp and ([TI|), Maxwell's equation is 

d^P P\\{t) 



u;{p,t)S{p,t) . 



(21) 



Ein{t) = -An(t) = -2e 



(2vr)3,;(P) 



f{P,t) + -v{P,t) 



(22) 



It is important to observe that this form for the internal field has been employed exten- 
sively in the study of back-reactions. However, our contribution is to employ it in conjunction 
with a time- dependent external field, which allows for the exploration of a richer variety 
of phenomena. 



B. Renormalisation 



The boson and fermion currents are 

d^P P||(t) 



jinit) = eg± 



(2vr)3u;(P) 



(23) 



where the integrand depends on the solution of the kinetic equation, Eqs. (|T2D-(|T4D, which 
must be such as to ensure the integral is finite. Simple power counting indicates that 
admissible solutions must satisfy 



|P|->oo 1 

v{P,t), fiP,t) < 



(24) 



To fully characterise the asymptotic behaviour we employ a separable Ansatz 



k=0 



k=0 



P 



Substituting these in Eqs. (p!^)-(|T^) and comparing coefficients, using Py 
which are valid at large |P|, we find the leading terms 

/4 = ^.e'^E^it) , = \eE{t) , Z2 = \eE{t) , 



UJ{P) 



(25) 



(26) 



with all the lower-order coefficients being zero. Substituting these results in Eq. (|20| ) it is 
clear that the conduction current is convergent. However, there is a logarithmic divergence 
in the polarisation current, Eqs. (|I9|) and (pT]), which is apparent in Eq. (p3D, but that is just 
the usual short-distance divergence associated with charge renormalisation. We regularise 
the polarisation current by writing v = {y — v^P\\/uj^) + vzP\\/uj^, so that 
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(27) 



-9±e 
where 



d^P P||(t) 



f{P,t) + Uv{P,t) 



eE{t)P\\{t)' 
^uj\P) 



-e'E^{t)I^{A) 



9± f d^P Pnit) ( 



A^oo g± 



871^ 



In 



(28) 



with A a cutoff on \P\, which effects a regularisation equivalent to Pauh-Villars. Intro- 
ducing the renormahsed charge, fields and current: 



(29) 



with Z = 1/(1 + e2/±(A)), and noting that eE^{t) = enS^it) and eA^{t) = eRA^{t), 
Eq. (^71) becomes 



d^P P||(t) 



(30) 



-9±eR 



(2vr)3^(P) 



eflg^(t)P||(t) 
4cj4(P) 



This defines a properly renormalised equation for the fields. Our procedure is technically 
different from that employed elsewhere P,p!8| but yields an equivalent result. Subsequently 
all fields and charges are to be understood as renormalised. 



C. Numerical results 



Equations ([T^ ) - ([T^) together with Maxwell's equation, Eq. (pTj) , form a coupled system 



of differential equations. To solve it we first evaluate the internal current from Eq. (jn\) 



at the primary time-slice using the initial conditions for the distribution function. That, 
via Eqs. (|T5p and (|TBp, provides an electric field, which we use to calculate the momentum 
distribution from Eqs. (|12[)-(|14[). This procedure is repeated as we advance over our time- 
grid. We use a momentum grid with 200 transverse- and 400 longitudinal-points, a time-step 
dt = 0.005, and A = 50 in Eq. (pH]). All dimensioned quantities are expressed in units of m, 
the parton mass. 

Spontaneous particle creation occurs in the presence of a strong field under whose influ- 
ence the vacuum becomes unstable and decays. Herein we induce this by a time- dependent 
external field and compare three different field configurations. The fields vanish at t —oo, 
and at if: = to the magnitude of the field increases and eventually leads to particle creation. 
Configuration (i): For comparison with Refs. [P,P|,p!l[, we solve the set of equations as an 
initial value problem without an external field, using an initial value of the electric field that 
is large enough to cause pair production. Configuration (ii): We employ 

Aex = -Aob^ [t/b + ln(2 cosh(t/6))], Ee^t) = AqB [tanh(t/6) + 1] , (31) 

which is an electric field that "switches-on" at t ~ — 6 and evolves to a constant value, 2Aob, 
in an interval t ~ 2/6. Configuration (iii): Is an impulse field configuration: 
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= Ao[tanh(t/6) + 1], E,,(t) = -Ao[b cosh\t/b)]-\ (32) 

which is an electric field that "switches-on" at t ~ —2b and off at t ~ 2b, with a maximum 
magnitude of A^/b a.t t = 0. Once this field has vanished only the induced internal field 
remains to create particles and affect their motion. 

For configuration (i) we fix an initial value of E{t = 0) = 10, in units of m, with 
= 4, and for bosons obtain the electric field and current depicted in Fig. |l|, where plasma 
oscillations are evident. The frequency of these oscillations increases with the magnitude of 
the field. The current exhibits a plateau for small t, when the particles reach their maximum 
velocity. This current opposes the field, and leads to a suppression of particle production 
and a deceleration of the existing particles. The effect of this is to overwhelm the field and 
change its sign with a consequent change in the direction of the particles' collective motion. 
The process repeats itself, yielding the subsequent oscillations that persist in the absence of 
additional interactions, such as collisions or radiation. The structure visible at the peaks and 
troughs of the current is not a numerical artefact. It is related to the field-strength/mass 
ratio, being more pronounced for large values, and occurs on a time-scale ~ Tqu, the vacuum 
tunnelling time, and hence can be characterised as Zitterbewegung. It disappears if an ideal- 
Markovian approximation to the source term is used because that cannot follow oscillations 



on such small time-scales [Q. The t-dependence of the p = distribution function is 
depicted in Fig. 0, where the beat-like pattern is the result of back-reactions and the rapid 
fiuctuations coincide with the Zitterbewegung identified in the current. 

Configurations (ii) and (iii) are alike in that the field "switches-on" at a given time. 
However, for (ii) the external field remains constant as t increases whereas in (iii) it "switches- 
off" after t ~ 2/b. The electric field and current obtained for bosons in these cases are 
depicted in Figs. |^ and ^: plasma oscillations are again evident. We plot the total electric 
field and thus it is evident in Fig. |^ that the internal electric field evolves to completely 
compensate for the persistent external field, which alone would appear as a straight-line at 
E(t) = 7. Unsurprisingly, as we see from Fig. H, a stable state is reached more quickly in 
the absence of a persistent electric field. Outside the temporal domain on which the vector 
potential acts, the initial value and impulse solutions are equivalent. 

We illustrate the results for fermions in Fig. ^ using the impulse configuration. The 
amplitude and frequency of the plasma oscillations are significantly larger than for bosons 
in a configuration of equal strength. Further, the stable state is reached more quickly 
because Pauli blocking inhibits particle production; i.e., no particles can be produced once 
all available momentum states are occupied. Pauli blocking also guarantees /_(p, t) < 1, for 
all t. 

We have also calculated the P\\- and p_L-dependence of / for both bosons and fermions. 
We find /+(t = 0) = 0; i.e., bosons cannot be produced with zero kinetic momentum, an 
effect readily anticipated from Eq. (^). For small t, f±{p,t) is a slowly varying function of 
p on its domain of support. However, with increasing t, f±{p,t) develops large-magnitude 
fluctuations without increasing that domain. The momentum- space position of the midpoint 
of the domain of support oscillates with a t-dependence given by the kinetic momentum: 
P|l=P||-eA(t). 

One additional observation is important here. The magnitude, oc Aq, of the electric 
flelds we have considered is large and hence the time between pair production events, Tprod, 
is small, being inversely proportional to the time-average of the source term, 5*. The period 
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of the plasma oscillations, Tpi, also decreases with increasing Aq but nevertheless we always 



have Tprod ^ Tpi- Thus, in contrast to the effect it has on the production process [p!7| , 
the temporal nonlocality of the non-Markovian source term is unimportant to the collective 
plasma oscillation. 



D. Collisions 

In the previous subsection we ignored the effect of collisions when treating the spon- 
taneous production of charged particles and subsequent evolution of the plasma. Now we 
consider the effect of a simple collision term 



c^^p^^t^ = i±m^i±m ^ (33) 

Tr 

where is the "relaxation time" and are the thermal equilibrium distribution functions 
for bosons and fermions: 

f±\p, t) = r f-,l/Tm^^^ ■ (34) 

exp[uj{p,t)/T{t)\ T 1 

Here T(t) is the "instantaneous temperature", which is a model-dependent concept p|, p!9| , p0| , 
and since our results are not particularly sensitive to details of its form we employ a simple 
parametrisation 

T{t) = T,, + (T„ - Te,) e-*'/*o , (35) 

with an equilibrium temperature T^q = 1.0, a maximum temperature Tm = 2.0, and a 
profile- width = 10 ~ Tpi. The collision term is added to the r.h.s. of Eq.(0), which 
becomes 

^^^ = S^ip,t) + C^ip,t). (36) 

This "relaxation time" approximation assumes that the system evolves rapidly towards 
thermal equilibrium after the particles are produced. It has been used before, both in the 



absence of back-reactions and including them [0,^,^, but with source terms that 
neglect fluctuations on short time-scales. 

We note that in the low density limit: f{p,t) <C 1, one can neglect the distribution 
function in the source term, Eq. (|), and Eq. ( ^6l) has the simple solution 



fim = f dt'exp (Si{p,t') + kSPlll) . (37) 

J —oo L Tr \ \ Tr J 

In our numerical studies we treat parameter and study the effect of C on the 

plasma oscillations. Our results for bosons using this crude approximation are depicted in 
Fig. Ij. For Tr ^ Tpi the oscillations are unaffected, as anticipated if l/r^ is interpreted as a 
collision frequency. For Tr ~ Tpi the collision term has a significant impact, with both the 
amplitude and frequency of the plasma oscillations being damped. There is a below which 
no oscillations arise and the systems evolves quickly and directly to thermal equilibrium. 
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IV. SUMMARY AND CONCLUSION 



We have studied spontaneous particle creation in the presence of back-reactions and 
colhsions, both of which dramatically affect the solution of the kinetic equation. The back- 
reactions lead to plasma oscillations that are damped by the thermalising collisions if the 
collision frequency is comparable to the plasma frequency. In electric fields where the period 
of the plasma oscillations is large compared to the time-scales characterising particle pro- 
duction, the non-Markovian features of the source term play little role in the back-reaction 
process. 

Plasma oscillations are a necessary feature of all studies such as ours but are they relevant 
to the creation of a quark-gluon plasma? If we set the scale in our calculations by assuming 
that fermions are created in the impulse configuration with (ej,) ~ Aqcd ~ 0.5^/a, the 
QCD string tension, then Aq = 10 with = 5 corresponds to an initial field strength 
eE ~ 15 (J and energy density ^E^ ~ 20^2. These are very large values but even so the 
plasma oscillation period is still large: Tpi = 5fm/c, and collisions can only act to increase 
that. We therefore expect that a quark-gluon plasma will have formed and decayed well 
before plasma oscillations can arise. On these short time-scales non-Markovian effects will 
be important. 

Our estimate shows back-reactions to be unimportant on small time-scales but that is 
not true of collisions. However, it is clear that in QCD applications they must be described 
by something more sophisticated than the "relaxation-time" approximation. 

Finally, with 1/Aqcd setting the natural scale, the finite interaction volume will clearly 
be important and hence the assumption of a spatially homogeneous background field must 
also be improved before calculations such as ours are relevant to a quark-gluon plasma. 
Ref. [p2| is a step in that direction. 
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FIG. 1. Time evolution for bosons of the total electric field and total current with initial values 
E = 10 and = 4. (Here, as in the text, all dimensioned quantities are given in units of the 
mass-scale m.) 
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FIG. 2. Time evolution of the distribution function /(O, t) with initial values E = 10 and = 4 
for bosons. The rapid fluctuations occur on the Vctcuuni tunnelling time-sccile, t^i^, and coincide 
with the Zitterbewegung identified in the current. In the absence of back-reactions these recurrent 



fluctuation packets are absent [17|. 
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FIG. 3. Time evolution for bosons of the electric field and the current for the step function 
external field, Eq. (^T|), with Aq = 14.0, b = 0.25 and the coupling = 6. 
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FIG. 4. Time evolution for bosons of the electric field and current for an impulse external field, 
Eq. (H), with ^0 = 10.0, b = 0.5 and the coupling = 5. 
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FIG. 5. Time evolution for fermions of the electric field and current for an impulse external 
field, Eq. (H), with = 10.0, b = 0.5 and the coupling = 4. 
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FIG. 6. Time evolution for bosons of the electric field obtained using different relaxation times 
the collision term of Eq. (|33|), and with the impulse external field, Eq. (^2|), where Aq = 7.0, 
= 0.5 and the coupling = 4. 
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